A proper description of clumping in hot star winds: 
the key to obtaining reliable mass-loss rates? 
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Abstract: Small-scale inhomogeneities, or 'clumping', in the winds of hot, massive stars are conventionally 
■ included in spectral analyses by assuming optically thin clumps. To reconcile investigations of different diag- 
nostics using this microclumping technique, very low mass-loss rates must be invoked for O stars. 

Recently it has been suggested that by using the microclumping approximation one may actually drastically 
underestimate the mass-loss rates. Here we demonstrate this, present a new, improved description of clumpy 
winds, and show how corresponding models, in a combined UV and optical analysis, can alleviate discrepancies 
between previously derived rates and those predicted by the line-driven wind theory. Furthermore, we show 
that the structures obtained in time-dependent, radiation-hydrodynamic simulations of the intrinsic line-driven 
' instability of such winds, which are the basis to our current understanding of clumping, in their present-day 
form seem unable to provide a fully self-consistent, simultaneous fit to both UV and optical lines. The reasons 
• for this are discussed. 

d 

1 Introduction 

The winds from hot, massive stars are described by the radiative line-driven wind theory (Castor et 
^ ■ al. 1975), in which the standard model assumes the wind to be stationary, spherically symmetric, and 
homogeneous. Despite its apparent success (e.g., Vink et al. 2000), this model is probably oversim- 
plified. In particular, much evidence for a time-dependent, small-scale inhomogeneous wind (that is, 
a 'clumped' wind) has over the past years accumulated, from the theoretical as well as the observa- 
tional side (for an overview, see Hamann et al. 2008). In the following, we investigate and discuss the 
indirect evidence for wind clumping that has arisen from quantitative spectroscopy aiming to infer 
mass-loss rates from observations. 



2 Deriving empirical mass-loss rates from hot star winds 

The main mass-loss diagnostics of OB-star winds are UV resonance lines, H a line emission (and 
other recombination lines), and infra-red and radio continuum emission. Recently, X-ray emission 
lines have also been added to the set (e.g., Cohen et al. 2010). Here we will focus on resonance 
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lines and H Q , discussing the influence of optically thick clumping on these diagnostics and using the 
well-studied Galactic O supergiant A Cep as a test bed. 



2.1 Smooth wind models and microclumping 

When smooth wind models are used, mass-loss rates inferred for a given star, but from different 
diagnostics, can vary substantially. For example, in an analysis by means of the unified NLTE model 
atmosphere code fastwind (Puis et al. 2005), we find that the phosphorus v (Pv) UV resonance 
lines suggest a mass-loss rate approximately 20 times lower than the one required for a decent fit 
of the H Q emission. That is, depending on which diagnostic is used, the inferred mass-loss rate of 
A Cep can vary by more than an order of magnitude. This inconsistency has been interpreted as a 
consequence of neglecting clumping when deriving these rates. 

Wind clumping has traditionally been included in diagnostic tools by assuming statistically dis- 
tributed optically thin clumps and a void inter-clump medium, while keeping a smooth velocity field. 
The main result of this microclumping approach is that mass-loss rates derived from smooth models 
and diagnostics that depend on the square of the density (such as H Q in OB-star winds) must be scaled 
down by the square root of the clumping factor, / c i(r) = (p 2 )/ (p) 2 with the angle brackets denoting 
spatially averaged quantities. On the other hand, processes that depend linearly on the density (such 
as the UV resonance lines) are not directly affected by microclumping. Thus, with such an optically 
thin clump model, the only way to reconcile the above large deviations in rates would be to accept 
the very low mass-loss rate indicated by PV, and so assume very high clumping factors (in this case 
f cl m 400!) to simultaneously reproduce H Q . Using this technique, we would derive a mass-loss rate 
for A Cep that is roughly an order of magnitude lower than predicted by the line-driven wind theory 
(according to the mass-loss recipe in Vink et al. 2000). Similar results have been found by, e.g., 
Bouret et al. (2005) and Fullerton et al. (2006). 

Such low mass-loss rates would have enormous implications for massive star evolution and feed- 
back, and of course also cast severe doubts on the validity of the line-driven wind theory. In addition, 
the extreme clumping factors we had to invoke to reproduce H a are in stark contrast with theoretical 
predictions (e.g., Runacres & Owocki 2002, see also Fig. 0). A possible solution of this dilemma 
is that clumps are not optically thin for the diagnostic lines under consideration, which would lead 
to underestimated rates if microclumping were still assumed (Oskinova et al. 2007; Owocki 2008; 
Sundqvist et al. 2010a). 



2.2 Relaxing the microclumping approximation 

Clump optical depths. The critical parameter determining the validity of the microclumping ap- 
proximation is the clump optical depth (which of course is different for different diagnostics). Ac- 
tually, for both resonance lines and H a , we may quite easily obtain reasonable estimates for this 
quantity. In the following, all radii are given in units of the stellar radius and all velocities in units 
of the terminal speed. The radial clump optical depth for a resonance line may then, in the Sobolev 
approximation, be written as 
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where q is the ionization fraction of the considered ion and / ve i a velocity filling factor (Owocki 
2008), defined as the ratio between the clump velocity span, 8v, to their velocity separation, in full 
analogy with the traditional volume filling factor. / ve i largely controls how a perturbed velocity field 
affects the line formation, and may be approximated by / ve i ps {Sv/Svplf^ 1 (Sundqvist et al. 2010b), 
where 5v@ is the clump velocity span of a corresponding model with a smooth velocity field. As 



seen from the middle term in Eq. [T] f ve \ also controls how the clump optical depth differs from the 
corresponding smooth one. n is a line-strength parameter proportional to the mass-loss rate and the 
abundance of the considered element, and taken to be constant within the wind (e.g., Puis et al. 2008). 
As an example, assuming a mass-loss rate 3.0 x 10~ 6 M /yr and a solar phosphorus abundance, we 
get k m 3.0 for the blue component of Pv in A Cep. Note that Eq.[T]assumes that the clumps cover a 
complete resonance zone. This is reasonable for the wind's inner parts, but may be questionable for 
its outer, slowly accelerating parts (see Sundqvist et al. 2010b). For our discussion here, however, a 
Sobolev treatment suffices. For H Q , the analogy to k is the parameter A (Puis et al. 1996^ which 
is proportional to the mass-loss rate squared and to the departure coefficient of the lower transition 
level. We may write the radial Sobolev clump optical depth for H a as 
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Assuming unity departure coefficients, and a mass-loss rate as above, A fa 1.5 x 1CT 3 for H Q in A Cep. 
Note the extra f c \ term in this expression, stemming from the p 2 -dependence of this diagnostic. In 
Fig [U we plot the H Q and Pv clump optical depths, using k and A as just given and f c \ = 9, which 
is a typical value. For simplicity, we also used q = 1 for the Pv line formation and a smooth velocity 
field, v = 1 — 0.99/r (i.e., a '/3=V law), which gives f c \=fvei~ 1 , when calculating the clump optical 
depths displayed in the figure. Clearly, within these approximations the clumps remain optically thick 
for Pv throughout the entire wind. For H Q , the p 2 -dependence makes the optical depth decrease faster 
with increasing radii, so this line is optically thick only in the lower wind. Based only on these simple 
estimates, we may therefore expect that resonance lines such as Pv should be sensitive to deviations 
from microclumping in the complete line profile, whereas recombination lines such as H Q should be 
affected only in the line core. We note also that even if we were to reduce the ion fraction of Pv to, 
say, q = 0.1, the clumps would still be optically thick. This illustrates the necessity of relaxing the 
microclumping approximation for typical mass-loss line diagnostics of hot star winds. 




Figure 1 : Clump optical depths for H Q and Pv line formation as functions of wind velocity. Stellar 
and wind parameters as for A Cep, see text. 



Wind models. To investigate effects on the line formation from optically thick clumps, a non- 
void inter-clump medium, and a non-monotonic velocity field, we create (pseudo-)2D and 3D wind 
models, following the 'patch method' of Dessart & Owocki 2002. In this technique, the full wind 
is 'patched' together by assembling inhomogeneous, spherically symmetric, wind-snapshots in ra- 
dially independent slices. We construct winds using both self-consistent, time-dependent, radiation- 
hydrodynamic (RH) simulations (computed following Feldmeier et al. 1997) and empirical, stochastic 
models. Line synthesis of resonance lines is carried out using the Monte-Carlo method described in 
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Sundqvist et al. (2010a), whereas for H Q we have developed a new radiative transfer code, in which 
we solve the 'formal integral' within our 3D winds, assuming that the departure coefficients are un- 
affected by optically thick clumping, which should be reasonable for the O-star winds discussed here 
(Sundqvist et al. 2010b). Hydrogen departure coefficients, and ionization fractions of Pv, are calcu- 
lated with fastwind, accounting for radially dependent microclumping. Stellar rotation is treated 
by the standard convolution procedure of a constant v smi (thus neglecting differential rotation), set 
to 220 km s _1 . 

When creating our stochastic models, we take an heuristic approach and use a set of parameters 
to define the structured medium. If clumps are optically thick for the investigated diagnostic, the 
line formation will generally depend on more structure parameters than just f c \. These parameters 
(for a two component medium) were defined and discussed in Sundqvist et al. (2010a). They are the 
inter-clump medium density, the physical distances between the clumps, vp5t, set by the time interval 
St between the release of two clumps and in our geometry equal to the porosity length (Owocki et 
al. 2004), and finally the velocity filling factor, / vel (see previous paragraph). We stress that these 
parameters are essential for the radiative transfer in an inhomogeneous medium with optically thick 
clumps, and not merely 'ad-hoc parameters' used in a fitting procedure. We notice also that in our 
stochastic models they are used to define the structured wind, and so are independent of the origin to 
the inhomogeneities. This should be distinguished from the RH simulations, in which the structure 
arises naturally from following the time evolution of the wind and stems directly from the line-driven 
instability. For these simulations then, the averaged structure parameters are an outcome. 



3 Results from an exemplary study of A Cep 
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Figure 2: Left and middle panels: Observed and modeled Pv and H a line profiles in A Cep. Solid 
black lines are calculated from empirical, stochastic wind models, and blue dashed ones from RH 
simulations. The dashed-dotted line in the left panel illustrates corresponding results using the micro- 
clumping approximation. Observations (red dotted lines) are from Fullerton, et al. (2006) (Pv) and 
Markova et al. (2005) (H a ). Right panel: Corresponding clumping factors. 



We have carried out a combined Pv and H a study, using both RH models and stochastic ones. 
We find that synthetic spectra computed directly from the RH models are unable to reproduce the 
diagnostic lines (Fig. 0). Two main problems are identified: i) the absorption toward the blue edge 
of Pv is too deep, and ii) the core emission of H Q is much lower than observed. Note that changing 
the mass-loss rate for which the RH model is calculated (M=1.5 x lO~ 6 M /yr) does not resolve 
these issues, or even alleviate them; if for example a higher rate is adopted, the wings of H a become 
much too strong. Consequently, we apply our stochastic models, aiming to empirically capture the 
essence of the structured wind. By means of these models, we obtain consistent fits essentially by 
increasing the clumping in the lower wind, but also by adopting somewhat lower velocity spans 



of the clumps. The second point largely resolves the issue of reproducing the observed Pv lines, 
and was extensively discussed in Sundqvist et al. (2010a). Regarding the first point, the observed 
H a absorption trough followed by the steep incline to rather strong emission can only be reproduced 
by our models if clumping is assumed to start at a velocity only marginally lower than predicted by 
the RH models (see also Puis et al. 2006; Bouret et al. 2008), but with a much steeper increase with 
velocity (Fig.|2l right panel). Regarding the outermost wind, RH simulations by Runacres & Owocki 
(2002), which extend to much larger radii than those used here, indicate that the clumping factor 
there settles at approximately four. f d ^4 is consistent with our derived mass-loss rate (see below) 
and the constraints from radio emission derived by Puis et al. (2006), suggesting that the outermost 
wind is better simulated than the inner by current RH models. However, let us point out that the use of 
the so-called smooth source function formalism in our RH models probably leads to overly damped 
perturbations in the inner wind (e.g., Owocki & Puis 1999), which might at least partly explain the 
discrepancies between theoretical clumping factors and those inferred from observations. 

The basic expectations from the clump optical depth estimates in Sect. 12.21 are confirmed by our 
detailed analysis. The Pv line profiles are clearly weaker than those calculated using microclumping 
(Fig.[2l left panel), whereas H a is affected only in the line core (not shown). Finally, our derived mass- 
loss rate for A Cep is the same as the one used in our RH simulations. It is approximately two times 
lower than predicted by the line-driven wind theory, but a factor of five higher than the corresponding 
rate derived assuming microclumping. This suggests that only moderate reductions of current mass- 
loss predictions for OB-stars might be necessary, and illustrates how a correct description of clumping 
is pivotal for obtaining consistent, reliable estimates of mass-loss rates. 

Acknowledgements 

J.O.S and J.P acknowledge a travel grant from the DFG cluster of excellence. 

References 

Bouret, J.-C, Lanz, T., & Hillier, D. J. 2005, A&A, 438, 301 

Bouret, J., Lanz, T., Hillier, D. J., et al. 2008, in Clumping in Hot-Star Winds, ed. W.-R. Hamann, A. Feldmeier, 

&L. M. Oskinova, 31 
Castor, J. I., Abbott, D. C, & Klein, R. I. 1975, ApJ, 195, 157 
Cohen, D. H., Leutenegger, M. A., Wollman, E. E., et al. 2010, MNRAS, 405, 2391 
Dessart, L. & Owocki, S. P. 2002, A&A, 383, 1113 
Feldmeier, A., Puis, J., & Pauldrach, A. W. A. 1997, A&A, 322, 878 
Fullerton, A. W., Massa, D. L., & Prinja, R. K. 2006, ApJ, 637, 1025 
Hamann, W.-R., Feldmeier, A., & Oskinova, L. M., eds. 2008, Clumping in hot-star winds 
Markova, N., Puis, J., Scuderi, S., & Markov, H. 2005, A&A, 440, 1 133 
Oskinova, L. M., Hamann, W.-R., & Feldmeier, A. 2007, A&A, 476, 1331 

Owocki, S. P. 2008, in Clumping in Hot-Star Winds, ed. W.-R. Hamann, A. Feldmeier, & L. M. Oskinova, 121— 

Owocki, S. P., Gayley, K. G., & Shaviv, N. J. 2004, ApJ, 616, 525 

Owocki, S. P. & Puis, J. 1999, ApJ, 510, 355 

Puis, J., Kudritzki, R.-P, Herrero, A., et al. 1996, A&A, 305, 171 

Puis, J., Markova, N„ Scuderi, S., et al. 2006, A&A, 454, 625 

Puis, J., Urbaneja, M. A., Venero, R., et al. 2005, A&A, 435, 669 

Puis, J., Vink, J. S., & Najarro, F. 2008, A&A Rev., 16, 209 

Runacres, M. C. & Owocki, S. P. 2002, A&A, 381, 1015 

Sundqvist, J. O., Puis, J., & Feldmeier, A. 2010a, A&A, 510, Al 1 

Sundqvist, J. O., Puis, J., Feldmeier, A., et al. 2010b, submitted to A&A 

Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2000, A&A, 362, 295 



